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Abstract 

We consider a scale invariant model which includes a term in action and show that a stable 
"emerging universe" scenario is possible. The model belongs to the general class of theories, where 
an integration measure independent of the metric is introduced. To implement scale invariance 
(S.I.), a dilaton field is introduced. The integration of the equations of motion associated with the 
new measure gives rise to the spontaneous symmetry breaking (S.S.B) of S.I. After S.S.B. of S.I. 
in the model with the term (and first order formalism applied), it is found that a non trivial 
potential for the dilaton is generated. The dynamics of the scalar field becomes non linear and 
these non linearities are instrumental in the stability of some of the emerging universe solutions, 
which exists for a parameter range of the theory. 
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I. INTRODUCTION 



In modern cosmology our notions concerning the early universe have introduced a new 
element, the inflationary phase of the early universe which provides an attractive sce- 
nario for solving some of the fundamental puzzles of the standard Big Bang model, like the 
horizon and the flatness problems as well as providing a framework for sensible calculations 
of primordial density perturbations. 

Even in the context of the inflationary scenario however one encounters the initial singu- 
larity problem which remains unsolved, showing that the universe necessarily had a begin- 
ning for generic inflationary cosmo logical model 12 1 . 

Recently models that can avoid those conclusions have been discovered pH9|. The way 
to escape the singularity in these models is to violate the geometrical assumptions of these 
theorems, which assume i) that the universe has open space sections ii) the Hubble expansion 
is always greater than zero in the past. In [3|,jJ] the open space section condition is violated 
since closed Robertson Walker universes with k = 1 are considered and the Hubble expansion 
can become zero, so that both i) and ii) are avoided. 

In jsl even models based on standard General Relativity, ordinary matter and minimally 
coupled scalar flelds were considered and can provide indeed a non singular (geodesically 
complete) inflationary universe, with a past eternal Einstein static Universe that eventually 
evolves into an inflationary Universe. 

Those most simple models suffer however from instabilities, associated with the instability 
of the Einstein static universe. The instability is possible to cure by going away from GR, 
considering non perturbative corrections to the Einstein's fleld equations in the context of 
the loop quantum gravityj^, a brane world cosmology with a time like extra dimension 6], 
considering the Starobinski model for radiative corrections (which cannot be derived from 
an effective action) exotic matterjSjf or f{R) theories in the presence of perfect fluids with 
w < 9|. In addition to this, the consideration of a Jordan Brans Dicke model also can 
provide a stable initial state for the emerging universe scenario jlol |. 

In this paper we propose a different theoretical framework where such emerging uni- 
verse scenario is realized in a natural way, where instabilities are avoided and a successful 
inflationary phase with a graceful exit can be achieved. 

We work in the context of a two measures theory (TMT) [ll| and more speciflcally in 
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the context of the scale invariant reahzation of such theories 12l4l5l|. These theories can 



provide a new approach to the cosmological constant problem and can be generalized to 



obtain also a theory with a dynamical spacetime [16| . We will show how the stated goals 
can be achieved in the framework of the TMT models. 

The paper will be organized as follows: In section II the effective field Eqs. in the 
Einstein frame, give rise to an effective potential for a dilaton field (needed to implement a 
model with global scale invariance) which presents a flat region. This model can describe 
an emerging universe scenario. In section III, we look at the generalization of this model 



15| by adding a curvature square, "i?^", and show that the resulting model contain a flat 
region. Furthermore the universe can evolve into an inflationary state which can undergo a 
successful graceful exit. We end with a discussion and conclusions section. 

II. THE EMERGING SCENARIO AFTER THE INTRODUCTION OF A 
TERM 

It is possible to obtain a model that through a spontaneous breaking of scale invariance 
can give us an emerging universe scenario. We consider the scale invariant action of the 



form 



11 



Sl = j Li<^d^x + j LiV^d^x + e j (g'"' R^,^iT)Y^/^d'^x. (1) 
Here, Li and L2 are given by 

Li = —RiT,g) + \g^''d,<pdA - ^(0), (2) 

and 

L2 = f/(0), (3) 

respectively. Also, $ is defined as $ = e^""^^ e ahcdd adv^hda^ cd d where i{>a (a = 1, 2, 3, 4) 
is a 4-scalar fields. The last term in the previous action is not only globally scale invariant, 
but also locally scale invariant, that is conformally invariant. 

Let us see which are the equations of motion. The variation of the action respect to the 
4-scalar fields, v^a, yields to 

K^M = 0, (4) 
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where = e^"'"'^eabcddu^bda^cdi3^d, and since det{A^) = we get d^Li = 0, for $ 7^ 0, 

and thus Li becomes 

Li = —R{T, g) + ]-g^^-d,cpd,cP -V = M, (5) 

where M is an integration constant. 

The variation of the action with respect to the metric (7'^^ gives 

i?^,(r)(— + 2eRV^) + $^0,^0,. -heR^ + U{4>))^g^, = 0. (6) 

It is interesting to notice that if we contract this equation with g^^ , the e terms do not 
contribute. Solving the scalar curvature from this and inserting in the other e - independent 
equation Li = M we get the solution for the ratio of the measures i.e. x = = M+v'((f)) ■ 

We see that under a conformal transformation given by 

V = ^^^'>9f^'^ = ix- 2KtR)g^^, (7) 

the new metric g^^, defines the known Einstein frame. Equations (E]) can now be expressed 
in this frame as 

\ ^ 

Rtiu - -Qf, uR = -^T'^t^ , (8) 

where 
with 

The constant A that appears in this latter expression could be obtained from the original 
lagrangian, where in it is added a term proportional to XQ^, where Q is the conformal factor. 
This renders the action to be 



Seff,A = I \/-Qd X 

where 



--R{g)+p{<P,R) 

K 



(11) 



R) = - Veff. (12) 



with Veff given by Eq. (|T0 



Also, Eq. (|5]) expressed in terms of g°'^ becomes -^-R(r, (?) + (x — 2KeR) \g^''d^(t)dy(l) — V = 
M, where we have written V in place of V{(f)). This allows us to solve for R and thus we get 

-KjV + M) + ir''d^(l)d,(l)x 

Note that when we insert Eq. ( |T3l) into Eq. ( |TOl) . it depends on the derivatives of the 
scalar field. It acts as a normal scalar field potential under the conditions of slow rolling or 
low gradients and in the case the scalar field is near the region M + V = 0. 

In the scale invariant case, where V and U are given by V{(j)) = fie°'^ and U{(j)) = f^e'^'^^, 



respectively [l2|, it is interesting to study the shape of K// as a function of in the case 
when change a little bit. In this case, Vg// can be regarded as a scalar field potential. 
Then from Eq. (fT3l) we get R = —k,{V + M), which inserted into Eq. (JTOl) gives, 

4(efi:2(/ie°<^ + M)2 + f^e^»*) ^ ^ ^ 

In Fig{T]we have plotted the effective potential, Ve/f as a function of the scalar field, 
for M = —1, e = —1, /i = 1/2, /2 = 1, A = 1/10 and n = 1. We should mention here that 
the choices of the values of the different parameters was done in virtue that they respect 
the conditions for a stable emerging universe. Notice that it shows a flat region of positive 
vacuum energy for large 0, a minimum obtained at zero without fine tuning. 

In the limit a0 — t- oo the effective scalar potential, K//, becomes Ve/f — )■ 4(^eK'^f'^+f2) ~^ "^^ 

Notice that in all the above discussion it is fundamental that M ^ 0. If M = the 
potential becomes just a fiat one, V^ff = ^^^2fi_^_f^-<j + A everywhere (not only at high values 
of a0). All the non trivial features necessary for a graceful exit, the other fiat region 
associated to the Planck scale and the minimum at zero if M < are all lost. The case 
M 7^ implies that we are considering a situation with S.S.B. of scale invariance. 

We now want to consider the detailed analysis of emerging universe solutions. We start 
considering the closed Friedmann-Robertson- Walker metric, 

ds"" = df - a{tf (^Y^Ti + + sm^edct)'')^ , (15) 

where a{t) is the scale factor, t represents the cosmic time and due to homogeneity and 
isotropy we take the scalar field to be a function of the cosmological time only, i.e. = 0(t). 

We will consider a scenario where the scalar field is moving in the extreme right region, 
i.e. for asymptotic values of (0 — t- oo). In this case, the expressions for the energy density 



p and pressure p are given by, 

A 



and 



p = + + + 50'), (16) 



p=|(l + /.V)0'-(C^ + i^0'), (17) 



where the constants A, B and C are given by 



A 

/2 + K-'eft 



. ,1,2 . (18) 

(19) 



(21) 



and 



/i /i 1 



..n , l,2/.. =ff-4-3^ + A, (22) 



respectively. 

It is interesting to notice that all terms proportional to 0^ behave like " radiation" , since 

= ^/^'^^ - = - 5)0' = j^'e^\ (23) 

and 

= ^Kh4>' + 50^ = (^6/.^ + B)4>' = ^-^Kh4>\ (24) 

so that p^i = ^ is satisfied. This is also a consistency check, since the e terms are associated 
to a conformally invariant term (in the first order formalism) and in the region — )■ oo 
that symmetry should hold, at least for those contributions. The e terms therefore do not 
contribute to the trace of the energy momentum tensor in the region ^ oo. 

The equations that determines such static universe a{t) — — constant, d = and 
d = give rise to a restriction on 0o- Since d = and it is proportional to p + 3p, we must 
require that p + 3p = 0, which leads to 

{B - AK'^e)^l - + C = 0. (25) 
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This equation has two roots, the first being given by 

and the second one is 
where 

A = /2 + 12/2Ae+12A/iV. 
From these expression we could get some specific range of the parameters in order that 

• 2 

either po and (f)Q become positive. In the following section we will restrict these parameters 
in order to describe an emerging universe. 



-2^2 + VA 
3//2e 



(26) 



(27) 



III. STABILITY OF THE STATIC SOLUTION IN THE SCALE INVARIANT 
MODEL 

We will now consider the perturbation equations. Considering small deviations of (j) from 
the static emerging solution value 0o and also considering the perturbations of the scale 
factor a, we obtain, from Eq. ( lT6l) . that 

6p = A^o6<p + 4(5 + ^KhA),pl6<p. (28) 

At the same time Sp can be obtained from the perturbation of the Friedmann equation 
3(^ + iJ^) = up, which gives under a perturbation respect to the background static solution 

6 

^(5a = nSp. (29) 

On the other hand, from the second order Friedmann equation, "^"'""^j"^"" = —up, we get 
that 

2 2 

— = -2kpo = -upo = ^loi^Po, (30) 

where we have used po = — po/3 at a = oq, and we have chosen to express our result in terms 
of Qq, defined by po = {Qq — l)po, which for the emerging solution has the value Qq = |. 
Using this in Eq. (l29l) . we obtain 

5p = ^^5a. (31) 

ao 
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Substituting this latter expression into Eq. ( I28ll yields to a linear relation between 50 and 
5a so that 

(50 = Do^a, (32) 

where 



ao<P^{A + A{B + \AK-'ecPl)) 
Since we could write that p = [Vt — l)p, with 



= 2 1 



(33) 



(34) 



where, 



V;/; = C + 50^ (35) 
and therefore, the perturbation of the second order Friedmann equation leads to, 

-'^ + 2— = -K5p = -K5{{n-l)p), (36) 

to evaluate this, we use ( IMI) . ( l35i) and the expressions that relate the variations in a and 
fl32|) . Defining the "small" variable /3 as a(t) = ao(l + /3) we obtain, 2/3(t) + iyo/3(t) = , 
where. 



12 B0, 



A + (j)l (45 + 2/t2eA) Po 



(37) 



with Po = ^(l + «:'e05)f + C + 50^. 

Stability of the static solution requires that Wq > 0. This static solution has to respect 
the conditions (pl > and po > 0. Only the second solution, i.e. 02 (see Eq. flTH) ). fulfills 
these requirements, providing that the parameters satisfy the following inequalities 



e < (38) 

-ef^ < /2 (39) 

< A < ^ (40) 

12e ( 1 + ^ 



Here, we have taken k = 1. 
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FIG. 1: The form of effective potential Veff{4') versus the scalar field (f). We have used M = — 1, 
e = -1, /i = 1/2, /2 = 1, A = 1/10 and k = 1. 

IV. STABILITY OF THE STATIC SOLUTION: DYNAMICAL SYSTEM. 

The study of the stability of the static solution and the properties of the different equilib- 
rium points could be done in a more systematic way by using a dynamical system approach. 
In this scheme we rewrite the Friedmann and the conservation of energy equations as an 
autonomous system in terms of the variables H and x = cf)^. In order to do so, we differen- 
tiating the Fridmann equation and after using the expressions for p and p, Eq. (16, 17), and 
the energy conservation equation we obtain: 

H =\--A{l + K'ex)x. (41) 

We can rewrite Eq. (14T]) and the energy conservation equation as the following two equa- 
tions. 
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H 



- C + Bx^ - A{1 + ex)x 
3 - 



(42) 



X 



3A(1 + K^ex)x 



H , 



(43) 



q + An^ex + 2Bx 



where we have ehminated a from Eq. fHTl) by using the Friedmann equation and we have 
used the definition of the variable x in order to get Eq. fH3l) from the energy conservation 
equation. The equations fH2!) and (l43l) are a two-dimensional autonomous system on the 
variables H and x. 

In order to study the stability of the static solutions we look for critical points of the 
system (^2]) and These points are 



The critical points have different properties depending on the values of the parameters of 
the model (/i, /2, e, A). At this moment we are not going to give an exhaustive description 
of these properties for all the critical points, instead, we are going to focus on the particular 
critical points which are related with static universe. From the definition of the variables H 
and X we can note that only the first two critical points Eqs. f l44[ correspond to a static 
universe. In order to study the nature of these two critical points we linearize the equations 
f l42|) and fH3|) near these critical points. From the study of the eigenvalues of the system we 
found that the first critical point, Eq. (144|) . could be a center or a saddle point, depending 
on the values of the parameters of the model. On the other hand the second critical point 
Eq. (H^ is a saddle. 

Stable static solutions correspond to a center. Then, by impose the critical point, Eq.( H^ . 
becomes a center, we recover the stability condition Eqs. (l38|) . ( 139|) and (j40|) . 




(44) 



(47) 



(45) 



(46) 
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In order to satisfy the requirements of stability we are going take the values e = — 1, 
/i = 1/2, /2 = 1 and A = 1/10. 

In Fig. [2] it is shown a phase portrait for four numerical solution to Eqs. f l42p and (1431) . 
Also, in this figure we have included the Direction Field of the system in order to have 
a picture of what a general solution look like. In this figure are four of the six critical 
points described in Eqs. (H^ - (H71) . One of this point is the center equilibrium point {H = 
0,x = 0.56), the saddle point {H = 0,x = 0.77), the point {H = 0.18, x = 1) is a future 
attractor and {H = — 0.18,x = 1) is a past attractor. The other equilibrium points are 
{H = 0.38, X = 0) and {H = —0.38, a; = 0) which are a future attractor point and a past 
attractor point respectively. 



H 




FIG. 2: Plot showing part of the Direction Field of the system and four numerical solutions. 



As we have mentioned, only critical points with H = represent a static universe. Then 
in our case only the critical point which is a center correspond to a solution which represent 
a static and stable universe. In Fig. [3] we show the Direction Field near this critical point 
{H = 0,x = 0.56) together with two numerical solution. 

V. INFLATION AND ITS GRACEFUL EXIT 

The emerging phase owes its existence to a strictly constant vacuum energy (which here 
is represented by the value of A) at very large values of the field (p. In fact, while for M = 
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FIG. 3: Plot showing the Direction Field near the center critical point and two numerical solution. 



the effective potential of tlie scalar field is perfectly fiat, for any M 7^ the effective potential 
acquires a non trivial shape. 

We consider then the relevant equations for the model in the slow roll regime, i.e. for 
small and when the scalar field is large, but finite and we consider the first corrections to 
the flatness to the effective potential. Dropping higher powers of <p in the contributions for 
the kinetic energy and in the scalar curvature R, we obtain 

P = ^7^' + 14//, (48) 

R= -k{V + M). (50) 

Here, as usual x = M+vi^) ■ the slow roll approximation, we can drop the second derivative 
term of and the second power of in the equation for and we get 

3i/70 = -v;'//> (51) 

3H^ = KVeff, (52) 

where V^jj = The relevant expression for V^jf will be that given by (fT4|) . i.e., where all 

higher derivatives are ignored in the potential, consistent with the slow roll approximation. 
We now display the relevant expressions for the region of very large, but not infinite 0, 
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these are: 



Veff = C + Ciexp(— a</)), (53) 

X = 2^exp{a<P) - 2M^, (54) 
Ji Ji 

and 

7 = 7o + 7iea;p(-a</)). (55) 

The relevant constants that will affect our results are, C, as given by fl22]) and Ci and 70 
given by 



and 



/2 + ^'eff 



respectively. 

Using Eq. ( 153|) we can calculates the key landmarks of the inflationary history: first, the 
value of the scalar field where inflation ends, (pend and a value for the scalar field 0=,, bigger 
than this (0=,, > (pend) and which happens earlier, which represents the "horizon crossing 
point". We must demand then that a typical number of e-foldings, like = 60, takes place 
between 0=,,, until the end of inflation at = (pend- 

To determine the end of inflation, we consider the quantity S = and consider the 
point in the evolution of the Universe where 5 = 1, only when 6 < 1, we have an accelerating 
Universe, so the point 6=1 represents indeed the end of inflation. Calculating the derivative 
with respect to cosmic time of the Hubble expansion using fl52|) and fl5T]) . we obtain that 
the condition 6 = 1 gives 

S = ^{V^ff/V.ffr = 1, (58) 

working to leading order, setting 7 = 70, Vejf = C and V^j^ = —aCiexpi^—acpend)^ this gives 
as a solution, 

exp{a(j)end) = -77^^=^, (59) 
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notice that if M and /i have different signs and if e < 0, Ci < for the allowed range of 
parameters the stable emerging solution, so — Ci represents the absolute value of Ci. We 
now consider 0* and the requirement that this precedes (pend by N e-foldings, 

N = j Hdt = j jd(j) =- J ^^^dcf), (60) 

where in the last step we have used the slow roll equation of motion for the scalar field f l5T|) 
to solve for 0. Solving H"^ in terms of Vg// using (152|) . working to leading order, setting 
7 = 7o and integrating, we obtain the relation between 0^, and 0end, 

exp{a(j)^) = exp{a(f)end) - ^ (61) 

Ck7o 

as we mentioned before Ci < for the allowed range of parameters the stable emerging 
solution, so that 0=,, > (pend as it should be for everything to make sense. Introducing Eq. 
fl59|) into Eq. fl6Tl) . we obtain, 

, , , Ci , a Na^ 

We finally calculate the power of the primordial scalar perturbations. If the scalar field 
had a canonically normalized kinetic term, the spectrum of the primordial perturbations 
will be given by the equation 

oc (63) 

I' 4- 

however, as we can see from ( HHj) . the kinetic term is not canonically normalized because of 
the factor 7 in that equation. 

In this point we will study the scalar and tensor perturbations for our model where the 
kinetic term is not canonically normalized. The general expression for the perturbed metric 
about the Friedmann- Robert son- Walker is 

ds^ = -(1 + 2F)dt^ + 2a{t)D^idx'dt + a^{t)[{l - 2ip)5ij + 2E^ij + 2hij]dx'dx\ 

where F, D, ip and E are the scalar type metric perturbations and hij characterizes the 
transverse-traceless tensor perturbation. The power spectrum of the curvature perturbation 
in the slow-roll approximation for a not- canonically kinetic term becomes Ref . [l7 1 (see also 



Refs.pjj) 



Mg^*.f;, (64) 
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where it was defined "speed of sound", Cs 5 cLS 



C 



2 



Px + 2XPxx' 

with P{X, 0) an function of the scalar field and of the kinetic term X = —{l/2)g'^'^d^(j)d,y(f). 
Here P x denote the derivative with respect X. In our case P{X, (p) = 7(0) X — Vg//, with 
X = 0^/2. Thus, from eq. (!64ll we get 

Ps = ki ^. (65) 

The scalar spectral index n<j, is defined by 

ns-l = = -26 - 77 - s, (66) 

a Ink 

where V = Th ^ ~ respectively. 

On the other hand, the generation of tensor perturbations during inflation would produce 

j 

gravitational wave. The amplitude of tensor perturbations was evaluated in Ref.pJ^], where 

2 /'2XPx-P' 

and the tensor spectral index n^, becomes 



f 2XP^x-P \ 

V ^Planck ) 



dink ' 
and they satisfy a generalized consistency relation 

Pt 

r = — = -8csnT. (67) 

Therefore, the scalar field (to leading order) that should appear in Eq. fl63l) should be 
^/7o0 and instead of Eq. (l65l) , we must use 

^ = J^. (68) 

The power spectrum of the perturbations goes, up to a factor of order one, which we will 
denote ki as {5p/ pY, so we have. 
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FIG. 4: The plot shows r versus for three values of a. For a = 1 solid line, a = 0.1 dash line 
and a = 0.01 dots line, respectively. Here, we have fixed the values M = —1, e = —1, /i = 1/2, 
/2 = 1, A = 1/10 and k = 1, respectively. The seven-year WMAP data places stronger limits on 



the tensor-scalar ratio (shown in red) than five-year data (blue) 



this quantity should be evaluated at = 0* given by ( l62i) . Solving for from the slow roll 
equation flSTj) . evaluating the derivative of the effective potential using fl53|) and solving for 
H from (!52|) . we obtain, to leading order, 



■exp(2a;0*), 



using then (1^ for exp(a0^,), we obtain our final result. 



Ps = h 



k'C , 1 



+ 



Na 



(70) 



(71) 



3 a/2 Vto^' 

it is very interesting first of all that Ci dependence has dropped out and with it all depen- 
dence on M. In fact this can be regarded as a non trivial consistency check of our estimates, 
since apart from its sign, the value M should not affect the results. This is due to the fact 
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that from a different value of M (altliough with the same sign), we can recover the original 
potential by performing a shift of the scalar field (f). 

In FigJH we show the dependence of the tensor-scalar ratio r on the spectral index n^. 
From left to right a = 1 (solid line), a = 0.1 (dash line) and a = 0.01 (dots line), respectively. 



From Ref. 



19| . two-dimensional marginalized constraints (68% and 95% confidence levels) 



on inflationary parameters r and n^, the spectral index of fluctuations, defined at ko = 0.002 
Mpc~^. The seven-year WMAP data places stronger limits on r (shown in red) than five- 
year data (blue)[20|. In order to write down values that relate Us and r, we used Eqs. (l66|) 
and f E7|) . Also we have used the values M = —1, e = —1, /i = 1/2, /2 = 1, A = 1/10 and 
K, = 1, respectively. 

From Eqs. flBU]) . (IMj) and flB7|) . we observed numerically that for a = 1, the curve r = r{ns) 
(see FigJl]) for WMAP 7-years enters the 95% confidence region where the ratio r ~ 0.011, 
which corresponds to the number of e-folds, ~ 32. For a = 0.1, r ^ 0.103 corresponds 
to A^ ^ 227 and for a = 0.01, r ~ 0.136 corresponds to A^ ^ 14137. From 68% confidence 
region for a = 1, r ~ 0.010, which corresponds to A^ ^ 34. For a = 0.1, r ~ 0.08 corresponds 
to A^ ~ 240 and for a = 0.01, r ~ 0.109 corresponds to A^ ~ 14279. We noted that the 
parameter a, which lies in the range 1 > a > 0, the model is well supported by the data as 
could be seen from FigJH 



VI. DISCUSSION AND CONCLUSIONS 



The consideration of a scale invariant action leads to exactly the type of scalar field 
potentials needed to obtain an emerging universe scenario. The addition of a i?^-term gives 
rise to solutions where e < 0, — e/f < /2 and < A < . a\ and to the stability of 

the initial Einstein universe solution. 

We find that this is accomplished because the eR^ term is shown to have the effect of 
introducing nonlinearities in the theory that stabilize the solutions. 

In addition in order for the whole idea to work we need to take M ^ 0, which means 
we are considering the S.S.B. of the theory and this is what gives the effective potential a 
non trivial shape, which allows a graceful exit from the inflationary phase and can give the 
adequate power of the density perturbations. 

It should be pointed out also that the theory studied here, in the context of an action 
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that contains a measure $ independent of the metric and with the use of the first order 
formahsm (the connections are considered as independent degrees of freedom in the action 
principle), leads to equations of motion that are only second order, i.e. only second deriva- 
tives of the metric and the matter fields appear (although higher powers of the derivatives 
of the scalar field do appear). This after the new measure $ and the connections are solved 
(through the equations of motion) in terms of the metric and matter fields and after we 
express our results in the Einstein frame. 

This is in contrast to the usual theories in the second order formalism and with 
standard measure everywhere in the action. In this case, these theories lead to fourth order 
equations for the metric field. In Ref. 2l| it was shown that the fourth order structure of the 
equations can be reformulated as a system of second order equations which contain however 
one additional degree of freedom, a scalar field. For studies of inflation in the usual B? 
theories in the second order formalism, see Ref. 22| . In the case of the usual B? theories this 



scalar field contains a potential with a fiat region, as shown in Fig{TJ 

Using the WMAP seven-year data, we have found some constraints for the parameters 
appearing in our model. In particular, FigJH shows that for the values of the parameter 
1 > a > 0, the model is well supported by the WMAP data. 
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